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LABORATOF OPERATIONS 

The  Laboratory  Operation*  of  The  Aerospace  Corporation  ia  conducting 
experimental  and  theoretical  investigation*  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver¬ 
satility  and  flexibility  have  been  developed  to  a  high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel¬ 
opments  is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 


Aerophysic*  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans¬ 
fer,  reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  laser*. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  specie*  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo¬ 
sensitive  materials  and  sensors,  high  precision  laser  ranging,  and  the  appli¬ 
cation  of  physic*  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics:  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi¬ 
conducting,  superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapon*  environment;  application  of  fracture  mechanics  to  stress  cor¬ 
rosion  and  fatigue-induced  fractures  in  structural  metals. 

Space  Physics  Laboratory:  Atmospheric  and  ionospheric  physics,  radia¬ 
tion  from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  In  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphero  the  effects  of  optical,  electromagnetic,  and  particulate  radia¬ 
tions  in  space  on  space  systems. 
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ABSTRACT 


The  chemical  and  energy  transfer  reactions  involved  in  the  helium- 
diluted  DF-CC>2  chemical  transfer  laser  system  are  discussed.  Experimental 
and  theoretical  rate  coefficient  literature  is  reviewed,  and  the  most  probable 
rates  are  selected  for  the  kinetic  model  of  the  reaction  system.  Results  of 
computer  experiments  show  the  relative  importance  of  these  reactions  for 
pulsed  laser  simulation.  Predicted  laser  performance  is  most  sensitive  to 
values  of  rate  coefficients  for  DF-CO2  transfer,  collisional  deactivation 
of  CO2(00°l)  by  DF,  and  collisional  deactivation  of  DF(v)  by  D  and  F.  A 
detailed  investigation  of  the  relationships  of  cavity  and  chemical  mecha¬ 
nisms  for  two  levels  of  initiation  is  presented  for  an  initial  composition  of 
XF:lF2:lD2:8CO2:40He  at  50  torr  and  300 *K  with  X  =  0.1  and  0.01.  For 
F  =  0.01,  a  260-psec  pulse  with  6%  chemical  efficiency  is  predicted,  while 
the  higher  level  initiation  level  (F  =  0. 1)  produces  a  52-psec  pulse  with  8% 
efficiency. 
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I.  INTRODUCTION 


Evaluation  of  the  performance  of  gas  and  chemical  lasers  requires  an 
understanding  of  the  kinetic  and  chemical  mechanisms  that  pump  and  deacti¬ 
vate  the  energy  levels  associated  with  lasing.  Normally,  a  simple  model 
may  be  constructed  that  explains  why  lasing  occurs  for  a  particular  species. 
However,  a  more  detailed  representation  of  the  mechanisms  is  needed  for  an 
accurate  prediction  of  laser  performance.  Recently,  comprehensive  theo¬ 
retical  models  for  computer  simulation  of  laser  performance  have  been 
developed  1 1-4  ]  that  utilize  a  detailed  formulation  of  the  kinetics  of  the 
lasing  system.  A  review  of  the  kinetics  for  the  CC^-^  laser  system 
is  given  by  Taylor  and  Bitterman  [5  ];  Cohen  [6  j  gives  a  similar  review 
for  the  H^+F^  chemical  laser.  With  the  discovery  and  development 
(7-16]  of  the  DF-CO^  chemical  transfer  laser,  it  is  desirable  to  assemble  a 
kinetic  model  that  will  facilitate  future  studies  of  this  laser. 

This  report  presents  the  results  of  our  efforts  to  gather  and  compile  a 
complete  set  of  the  required  rate  coefficients  based  on  data  available  in  the 
literature,  comparisons  with  similar  reactions  with  known  rate  coefficients, 
or  various  theoretical  approaches.  We  have  drawn  heavily  on  Refs.  [  5] 
and  (6];  we  include  data  prior  to  those  studies  only  when  necessary. 

A  schematic  of  the  DF-CO^  transfer  and  the  CO^  kinetics  is  shown  on 
an  energy  level  diagram  in  Fig.  1.  The  vibration -translation  (VT)  energy 
transfer  reactions  are  denoted  by  dotted  lines;  intramolecular  vibration- 
vibration  (VV)  reactions  are  shown  as  solid  lines;  and  the  intermolecular  VV 
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Vibrational  energy  level  diagram  for  CO?  kinetic*  including 
the  DF  transfer  energy.  The  states  (03^o)  and  (02^0)  are 
deleted  for  clarity.  These  state*  are  included  in  the  model. 


transfer  reactions  are  denoted  by  dashed  lines  with  arrows  to  the  final 

states.  Initial  states  of  these  intermolccular  VV  reactions  are  CO2(000(?) 

and  the  CO^  state  at  the  tail  of  the  dashed  arrow.  The  states  (03 *0)  and 

(03^0)  are  assumed  to  be  in  equilibrium  and  arc  statistically  combined  as 

2 

(030).  Reactions  involving  the  state  (02  0)  are  similar  to  those  for 
(02^0)  but  have  been  omitted  from  the  figure  to  preserve  clarity. 

We  limit  our  model  to  energy  levels  in  CO^  below  CC^fOO^l);  for  com¬ 
parison,  Fig.  1  shows  three  higher  states.  This  limitation  is  artificial; 
however,  we  expect  levels  above  (00®l)  to  play  a  minor  role  in  transfer  laser 
performance.  We  assume  DF  to  transfer  only  to  CC^  (00®1).  The  net  effect 
on  laser  performance  of  transfer  to  levels  higher  than  CC^fOO^l)  and  their 
ultimate  deactivation  to  (00®1)  should  be  small. 

The  laser  simulation  program  used  in  this  study  utilizes  rate  coeffi¬ 
cients  of  the  form  (2] 


k.  atd.-e/rt 

where  A  and  D  are  constants,  E  is  the  activation  energy,  R  is  the  gas  con¬ 
stant,  and  T  is  the  temperature.  All  rate  coefficients  are  presented  in  this 
form.  The  proposed  kinetic  model  is  used  with  this  program  to  determine 
the  relative  importance  of  kinetic  mechanisms  during  a  pulse  for  two  different 
initial  F  atom  concentrations.  This  model  has  Droven  successful  in  the  pre¬ 
diction  of  characteristics  of  the  HF  pulsed  laser  [17],  A  number  of  rate  coef¬ 
ficients  are  roughly  approximated;  therefore,  they  are  varied  to  determine 
the  sensitivity  of  laser  performance  to  their  values  under  typical  laboratory 
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conditions.  This  sensitivity  study  shows  that  the  peak  pulse  power  and  pulse 
energy  arc  sensitive  to  the  F^  +  D ^  pumping  reactions  and  to  the  temperature 
dependence  of  the  DF(v)-CC>2  transfer  rate  coefficient  and  its  variation  with 
v.  The  temperature  variation  of  the  rate  of  VT  deactivation  of  CC^OO®!)  by 
DF  and  the  DF-F  and  DF-D  VT  processes  are  also  important  in  predicting 
laser  performance.  For  the  reactant  mixtures  studied,  the  rate  of  equili¬ 
bration  of  the  lower  levels  in  CO^  is  so  fast  that  they  are  near  equilibrium, 
and  thus  modeling  is  insensitive  to  the  value  of  rate  coefficients  for  these 
levels. 

In  Section  II,  we  review  the  kinetics  for  the  D^  +  F^  system,  which  may 
be  used  for  modeling  the  DF  chemical  laser.  Sections  III  and  IV  review  the 
kinetics  of  DF  -CO^  vibrational  energy  transfer  and  the  kinetics  of  CC^ 
relaxation,  respectively.  Section  V  describes  the  computer  experiments, 
and  Section  VI  gives  a  brief  summary  of  the  work. 
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II.  THE  D2  +  F2  KINETICS 


Theoretical  and  experimental  studies  of  the  kinetics  of  H2  +  Cl2  and  of 
H2  +  F2  mixtures  have  been  made  to  test  proposed  theories  and  increase  our 
understanding  of  the  underlying  processes.  We  draw  heavily  on  these  and 
analogous  mixtures  for  modeling  of  D2  +  F2  kinetics  when  direct  experi¬ 
mental  measurements  are  not  available. 

A.  RECOMBINATION  -  DISSOCIATION  REACTIONS 

The  initiating  step  in  the  exothermic  chemical  reaction  is  usually 
dissociation  of  F2  or  D2>  Conversely,  the  recombination  of  D  and  F  with 
themselves  or  each  other  represents  a  potential  energy  loss  mechanism  for 
laser  systems.  Under  typical  operating  conditions,  these  reactions  have 
little  or  no  effect  on  laser  performance;  however,  we  include  them  for  com¬ 
pleteness.  For  the  dissociation  rate  of  F2,  we  use  the  rate  coefficient  given 
by  Ref.  [6 J, 


(1) 


F2+M^F  +  F+  M 


kj  =  5  x  1013e 


35.  30 


where  M  is  all  collision  partners.  All  rate  coefficients  are  given  in  units  of 
cc/mole-sec  and  0  =  -103/RT.  The  recombination  rate  data  of  Jacobs,  Giedt, 
and  Cohen  [18]  for  React.  (2)  are  used: 

(2)  D2(0)  +  M^D  +  D+  M 


-5- 


where 


k 


2 


10 


for  M  equals  all  species  except  and  D.  For  D ^  and  D,  this  rate  is 
increased  by  a  factor  of  1.  75  and  20,  respectively.  The  minus  subscript 
indicates  a  backward  rate.  The  populations  in  the  higher  levels  of  D£(v) 
are  small  and  are  not  included  in  React.  (2). 

The  rate  of  dissociation  of  HF  has  been  measured  by  Jacobs,  Giedt, 
and  Cohen  [19].  We  take  their  rate  coefficient  for  DF  dissociation  by  all 
species 

(3)  DF(v)  +  M  ^  D  +  F  +  M  v  =  0, . . . ,  9 

with 

_  1.2  x  1019  m-1.0_-(D*  -  EJ/RT 
k3  -  x  e  v 


where  D#  (137.  13  kcal)  is  the  dissociation  energy  of  DF(0),  Ey  is  the  vibra¬ 
tional  energy  above  ground  state,  and  n  is  the  number  of  vibrational  levels 
considered,  here  taken  as  10.  The  factor  i/n  is  used  so  that  the  total  recom¬ 
bination  rate  is  in  agreement  with  experiment;  thus,  all  rates  of  recombina¬ 
tion  have  equal  probability.  The  vibrational  energy  for  level  v  is  computed 
from  data  given  by  Herzberg  [20],  The  enthalpies  of  formation  are  listed  in 
Table  I  for  all  species  included  in  our  model. 


« 
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Table  I.  Enthalpy  of  Formation  of  Included  Species 


Species 

(kcal/mole) 

Species 

AH? 

(kcal/mole) 

DF(0) 

-64.98 

CO2(00°l) 

-87.331 

DF{1) 

-56.67 

CO2(ll10) 

-88.111 

DF(2) 

-48.62 

C02(030) 

-88.531 

DF(3) 

-40.83 

co2(io°o) 

-90.081 

DF(4) 

-33.29 

C02(0220) 

-90.236 

DF(5) 

-26.00 

C02(02°0) 

-90.381 

DF(6) 

-18.95 

CO  2(  01*0) 

-92. 141 

DF(7) 

-12.13 

C02(00°0) 

-94.051 

30 

Q 

-  5.56 

D,(0) 

fir 

DF(9) 

0.786 

2' 

d2(1) 

8.544 

F 

18.9 

D2(2) 

16.764 

13.  CHEMICAL  PUMPING  REACTIONS 


Thu  overall  rate  coefficients  of  the  chemical  pumping  reactions 

k4(v) 

(4)  F  +  d2  ^  DF(v)  +  D 
and 

k5(v) 

(5)  D  +  F2  DF(v)  +  F 

have  been  calculated  to  be  10^*  ^e and  10^*  by  Wilkins  [21  ] 

13  8  2  169 

using  absolute  reaction  rate  theory.  An  overall  rate  of  10  *  e  '  for 

has  been  obtained  by  Jaffe  and  Anderson  [22]  from  a  limited  number  of 
classical  trajectory  calculations.  In  addition,  Monte  Carlo  calculations  by 
Muckerman  [23]  indicate  that  the  overall  rate  for  React.  (4)  is  about  0.6 
times  that  of  F  +  H^,  in  agreement  with  Ref.  [22]  and  the  calculations  of 
Wilkins  [  24] . 

On  the  basis  of  heats  of  formation  in  Table  I,  these  reactions  release 
30.63  and  99.  33  kcal/mole,  respectively,  for  v  =  0.  The  distribution  of  this 
energy  among  the  vibrational  levels  of  DF(v)  for  React.  (4)  has  been  mea¬ 
sured  by  Parker  and  Pimentel  [25]  at  390*K,  using  an  equal  gain  chemical 
laser  experiment,  and  by  Schafer,  et  al.  [26]  and  Anlauf,  et  al.  [27],  using 
a  crossed  molecular  beam  and  a  discharge  flow  technique,  respectively. 
These  data  are  compared  in  Table  11.  The  data  of  Schafer,  et  al.  ,  are 
for  back  scattering  only  and  are  inconsistent  with  trends  found  for  the 
more  widely  studied  F  +  reaction.  Their  distribution  also  disagrees 
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Table  11.  F  +  Pumping  Distribution 


Vibrational 

E 

k 

k 

k 

V 

V 

V 

V 

Level 

(kcal/mole)a 

127] 

126] 

[25] 

0 

0 

---- 

0.012 

— 

1 

8.31 

0. 1 

0.019 

---- 

2 

16.36 

0.5 

0.051 

0.63 

3 

24.15 

1.0 

0.285 

1.0 

4 

31.69 

0.72 

1.0 

<1.19 

aVibrational  energy  above  ground  state 


with  the  laser  experiments  of  Basov,  et  al.  [14j.  The  pumping  distribution 
data  for  F  +  includes  the  equal  gain  measurement  at  539aK  of  Parker  and 
Pimentel  [25]’,  crossed  molecular  beam  studies  of  Schafer,  et  al.  [26J;  dis¬ 
charge  flow  system  measurements  of  Anlauf,  et  al.  127],  Polanyi  and  Tardy 
[28],  and  Jonathan,  et  al.  [29-31 J ;  and  the  trajectory  calculations  of  Wilkins 
(32).  These  data  are  compared  to  those  of  React.  (4)  in  Fig.  2.  The  abscissa 
in  the  figure  ((Ev  -  Eq)/|  AH||  is  the  ratio  of  the  vibrational  energy  of  a  given 
level  to  the  heat  released  for  the  particular  reaction.  From  the  figure,  we 
see  that  the  data  of  Anlauf,  et  al.  ,  are  the  most  consistent  with  the  expected 
trends;  we  use  this  distribution  and  the  total  rate  computed  by  Wilkins. 

Reaction  (4)  is  not  sufficiently  exothermic  to  produce  DF(v)  above 
v  =  4;  however,  there  is  some  evidence  in  the  system  [6]  that  this 

reaction,  running  in  reverse  at  the  higher  v  levels,  is  a  sink  for  excited 
species  produced  by  reactions  like  (5).  Hence,  we  use  rate  coefficients 
analogous  to  those  of  Cohen  for  Hg  +  F.,: 


k.4<v 


5)  =  4  X  10 


12-.0.  15  30 
T  e 


and 


k_4(v  *  6)  =  1.2  x  10I3T0,  15 


These  values  are  roughly  twice  as  large  as  predicted  by  the  recent  Monte 
Carlo  calculations  of  Wilkins  [33]. 

No  studies  have  been  made  to  determine  the  distribution  of  pumping 
by  React.  (5).  Evidence  that  this  reaction  pumps  to  v  =  9  is  given  by  the 


k(v)/k(v). 


1.5 


F  +  H2  -  OPEN  SYMBOLS 
F+D2  -  SOLID  SYMBOLS 


(Ev-E0)  /  I  AH  | 


Fig.  2.  Pumping  distributions  for  F  +  H2  and  F  +  D2:  equal  gain 
technique,  ORef.  [25];  crossed  molecular  beam  study, 
ORef.  [26];  discharge  flow  system,  V  Ref.  [27];  □  Ref. 
[28],  O  Refs.  [29]  and  f30j ,  O  Ref.  [31 J;  Monte  Carlo 
calculations,  D  Ref.  [32] . 


spectral  output  of  the  DF  laser  measurements  of  Basov,  et  al.  [14]. 

Jonathan,  et  al.  [29,  30)  have  reported  the  distribution  of  pumping  by  H  +  F^. 
We  expect  that  the  similarity  shown  in  Fig.  2  between  F  +  and  F  +  will 
hold  between  H  +  F^  and  D  +  F2.  In  Fig.  3,  we  plot  the  distributions 
reported  by  Jonathan,  et  al.  [29,  30).  By  connecting  the  data  points  with  a 
smooth  curve  and  denoting  the  vibrational  energy  levels  of  DF(v),  we  esti¬ 
mate  the  pumping  distribution  to  be 

k5(0):k5(I):k5(2):k5<3):k5(4):k5(5):k5(6):k5(7):k5(8):k5<9) 

=  0:0.16  :0.27  :0.40  :0.57  :0.72  :0.86  :  1  :0.97  :0.6l 

With  these  estimates  of  the  distribution  of  pumping  into  DF(v),  Reacts. 
(4)  and  (5)  convert  79%  and  46%,  respectively,  of  the  available  chemical 
energy  into  vibrational  energy. 

C.  VIBRATIONAL-TRANSLATIONAL  ENERGY  TRANSFER  IN  DF 

Vibrational  deactivation  of  an  excited  species  through  collisions  is  a 
significant  energy  loss  mechanism  in  chemical  lasers, 

(6)  DF(v)  +  M  ^  DF(v-l)  +  M 

This  rate  coefficient  has  been  measured  in  a  shock  tube  for  v  =  1  by  Bott  and 
Cohen  [34]  over  a  wide  temperature  range  for  M  =  H^,  N^,  Ar,  and  DF.  These 
measurements  of  the  Napier  time  pT  are  converted  to  rate  coefficients  by 
assuming  an  harmonic  oscillator  model  and  using  the  relation 
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(1) 


k  = 


RT 

pr,7?^ 


Stephens  and  Cool  [35 j  also  measured  the  DF-DF  VT  rate  using  a  chopped 
DF  laser  at  350° K.  These  experimental  rates  are  shown  in  Fig.  4  on  a 
Landau-Teller  plot  with  approximate  analytic  rate  coefficients  matched  to 
the  data  for  the  temperature  range  300*  to  2000*  K,  The  character  of  the 
DF-DF  rate  as  a  function  of  <  *mperature  is  similar  to  that  of  the  HF -HF 
rate  [6,  36 J.  The  rate  coefficients  matched  to  the  data  in  Fig.  4  are: 


6  =  3.  4  X  104T2,  2 

m  =  h2 

6a=  4.0  xloV'2 

M=  N2 

6b=  3.7X10-V-66 

M  =  Ar 

,  «  4.5  x  104T2-2  +  5.3  x  1016T'2,0 

M  =  DF 

Except  for  high-temperature  measurements  of  DF-F,  deactivation  of  DF(v)  by 
other  collision  partners  has  not  been  measured.  The  efficiency  of  in  deac¬ 
tivating  DF  is  close  to  the  efficiency  of  H £  in  deactivating  HF[34],  Hence, 
we  use  the  measurement  of  HF  relaxation  by  of  Bott  and  Cohen  [  36]  for 
relaxation  of  DF  by 

k6d  .  0.  14  T3-  66  M  =  D2 


We  estimate  the  efficiency  of  F2  in  deactivating  DF  to  be  approximately  the 
same  as  Ar,  and  the  rate  for  DF -He  is  taken  as  twice  the  DF -Ar  rate 
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2000  1000  600 


k =4.5x10  T  ♦  5.3xKTT 
DF-DF 


4  22 

k  =  3.4  x  10  T 


df-h2 


k  =  4  xio3 T2-2 

"DF-N2 

\  k  =  0.14  T366 

nhf-d2 


\  k s  3.7  x  I0"6T  4  66 


\° 

F-Ar 

so5 

_l _ 1 _ 1 _ 1 _ 1 _ 

-J _ 1 _ 1 _ 

0.06 

0.08  0.10  0.12  0.14  0.16 

•j*  ”1/3 

0.18  0.20  0.22  0.24 

Fig.  4. 

Measured  rate  coefficients  for 
0  Ref.  [34],  A  Ref.  [35]. 

VT  deactivation  of  DF(  1 ): 
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analogous  to  their  efficiency  in  deactiviating  HF  [36].  The  deactivation  rate 
for  OF  -  CO^  is  assumed  to  be  similar  to  that  of  For  all  these  reactions, 
we  assume  the  deactivation  rate  is  linear  with  v;  i.e.,  we  assume  an  har¬ 
monic  model.  There  has  been  speculation  regarding  the  variation  of  these 
rates  with  v;  however,  this  question  can  only  be  answered  by  careful 
experiment. 

We  expect  the  deactivation  rate  with  F  and  D  to  be  quite  fast.  Shock 

tube  measurements  of  the  HF  -F  VT  rate  near  2000*K  by  Bott  and  Cohen  [37] 

show  relaxation  by  F  to  be  26.  5  times  more  efficient  than  HF  self- relaxation 

at  2000°K.  Also  at  2000*K,  shock  tube  measurements  of  Blauer,  et  al.  [38] 

1 3 

show  that  HF-F/ DF -F  ~  2.4  and  DF  -  F  is  approximately  1.5  x  10  cc/mole-sec, 
which  is  in  agreement  with  the  shock  tube  measurements  of  Bott  and  Cohen  [34]. 
Ridley  and  Smith  [39]  have  measured  the  rate  coefficients  for  deactivation  of 
HCl(v)  by  Cl  for  v  =  3,  2,  and  1  to  be  4.8  X  1012,  1.9  X  1012,  and  5.8  X  1011 
cc/mole-sec  at  room  temperature.  We  have  used  these  rates  to  approximate 
the  efficiency  of  F  and  D  in  deactivating  DF(v).  Thus,  for  M  =  D  and  F,  we  take 


k6e<l> 


5.  0  x  10 


12  1.250 
e 


k6e<2> 


5.  0  x  10 


12  0.550 
e 


k6e(v>3)=  S.0X1012 

The  activation  energies  for  v  =  1  and  2  are  chosen  to  approximate  the  mea 
surements  of  Ref.  [39]. 
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These  values  are  consistent  with  extrapolation  of  Ref.  [38]  and  with 
recent  trajectory  calculations  of  the  HF  -H  VT  rate  by  Wilkins  [33]  where 
multiple  quantum  deexcitation  was  predicted.  We  have  deviated  from  the 
ky  v  j  =  vkj  q  form  for  these  reactions,  since  they  are  apparently  chemical 
in  nature. 

D.  VIBRATIONAL-VIBRATIONAL  ENERGY  TRANSFER  IN  DF 

Since  the  efficiency  of  quantum  exchange  for  Av  >  1  is  expected  to  be 
orders  of  magnitude  less  than  for  Av  =  1,  we  follow  Ref.  [6]  and  consider 
only 

(7)  DF(v)  +  DF(v')  ^  DF(v  + 1)  +  DF(v'-l) 

where  v  -  v'  =  0,  1,  2.  Cohen  includes  up  to  v  -  v1  =  3;  however,  with  the 
increased  complexity  of  this  system,  we  must  limit  the  number  of  reactions. 
At  present,  no  experimental  data  on  VV  transfer  in  DF  are  available.  For 
reactions  of  this  type,  the  rate  varies  with  degree  of  resonance  and,  hence, 
would  be  expected  to  vary  with  the  rotational  levels  of  the  colliding  molecules. 
We  make  no  attempt  to  correct  for  this  effect  and  choose  a  rate  suggested  by 
Cohen  [6]  that  is  based  on  the  theory  of  Rapp  and  Englander -Golden  [40]  and 
the  data  of  Chen  and  Moore  [41]  for  HC1 

kO,l;l,0=ZX4xl0'5T*,0'A‘'/433 

where  Z  is  the  gas  kinetic  bimolecular  collision  rate  and  A  v  is  the  transla¬ 
tional  energy  released  in  wave  numbers.  The  gas  kinetic  bimolecular 
collision  rate  is  given  by 
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(2) 


and  N^  is  Avogadro's  number.  Table  Ill  contains  a  list  of  the  collision 
diameters  oj  and  the  masses  mj  used  in  our  gas  kinetic  calculations.  We 
assume  that  the  rate  coefficients  behave  as  though  the  vibrational  levels  of 
the  molecule  are  harmonic  and 

kv,  V+  Uv',  v'-l  *  <v  +  1>V'Z  *  4  X  10'5t  *  10'Av/433 

where  Avis  the  translational  energy  released  for  each  reaction. 

Recent  measurements  [42]  for  HC1  VV  reactions  indicate  that  this 
theory  may  not  be  appropriate  for  hydrogen  halides  with  v  and  v'  >  1. 
Although  these  VV  rates  are  in  question,  they  probably  do  not  play  a  major 
role  in  predicting  laser  performance,  as  suggested  by  computer  calculations 
of  Kerber,  Emanuel,  and  Whittier  [43]  and  the  analytic  study  by  Emanuel  and 
Whittier  [44]  for  the  HF  pulsed  laser. 

E.  VIBRATIONAL-VIBRATIONAL  ENERGY  TRANSFER 
BETWEEN  DF  AND  D.. 

Since  vibrational  energy  transfer  between  DF  and  D2  may  represent  a 
significant  loss  of  laser  energy,  we  include  the  reaction 

(8)  DF(v)  +  D2(v')  DF(v-l)  +  D2(v'  +1) 

in  our  kinetic  model.  We  use  the  same  probability  for  the  backward  rate  of 
this  reaction  as  for  DF-DF  VV  reactions, 
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(3) 


k_g  =  (v*  +  l)vZ  x  4  x  10"5T  x  l£>-^/433 

Recent  measurements  by  Airey  and  Fried  C  45)  of  the  rate  of  the  HF-H^  VV 
reaction  are  of  the  same  order  as  predicted  by  Eq.  (3). 

F.  VIBRATIONAL- TRANSLATIONAL  ENERGY  TRANSFER  IN 

These  reactions  are  included  to  prevent  unrealistic  inversions  in  the 
D2  populations.  Hence,  we  include 

(9)  D2(v)  +  M  ^  D2(v-1)  +  M 

This  rate  is  determined  from  the  compilation  of  Napier  times  given  by 
Stretton  in  Ref.  [46].  For  M  =  Ar,  this  gives 

k  =  1.  51  x  10‘5T4,  6 

We  use  this  rate  for  all  species  except  D2>  D,  and  He.  For  M  =  D2,  we 
increase  the  rate  by  a  factor  of  4.  8  to  be  consistent  with  Ref.  [461,  and  we 
estimate  that  D  and  He  are  at  least  twice  as  effective  as  Ar. 

The  proposed  kinetic  model  for  the  F2  +  D2  chemical  system  is 
summarized  by  Reacts.  (1)  -  (9)  of  Table  IV. 
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Table  IV.  Recommended  Rate  Coefficients  for  the  D-,  -  F,  -  CO 
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III.  DF  -  C02  VIBRATIONAL  ENERGY  TRANSFER 


The  most  important  reaction  in  the  DF-C02  laser  system  is  the 
transfer  of  vibrational  energy  from  DF(v)  to  CO2(00^1) 


(10) 


DF(v)  +  C02(00°0)  ^  DF(v-l)  +  CO2(00°l) 


12 

The  rate  coefficient  of  this  reaction  has  been  determined  to  be  3.  82  x  10 
cc/mole-sec  for  v  =  1  by  Stephens  and  Cool  [3  5]  using  a  laser-induced  fluores¬ 
cence  technique  at  350*K.  Preliminary  measurements  of  Cool  [47]  at  300 *K 

-2 

suggest  that  this  rate  coefficient  may  vary  as  strongly  as  T  .  Basov,  et  al. 

[14]  estimate  this  rate  to  be  9  X  10**  cc/mole-sec  at  300*K  from  a  DF-C02 

pulsed  laser  experiment.  The  rate  of  vibrational  energy  transfer  from  other 

hydrogen  halides  to  C02  has  also  been  measured  by  the  laser  fluorescence 

technique.  Stephens  and  Cool  [35]  found  k^F =  8. 1  x  10**  cc/mole-sec 

2  12 

at  350°K,  and  Chen,  et  al.  [48]  report  k^ci-CO  =  *‘®  X  ***  cc/mole-sec 

12  2 

and  kpji.CX)  =  2.  4  x  10  cc/mole-sec.  Although  the  degree  of  resonance 
for  v  =  1  is  1612  cm"*,  537  cm'*,  and  -116  cm"*,  respectively,  for  these 
reactions  compared  to  558  cm’'  for  DF,  their  efficiencies  of  transfer  are 
all  the  same  order  of  magnitude  as  the  DF  -  C02  value. 

An  analysis  of  the  energy  spacing  of  the  Av  =  1  transitions  in  DF  (Fig.  1) 
shows  that  the  closest  resonance  for  J  =  0  is  at  v  =  8.  We  have  also  calcu¬ 
lated  the  relative  probabilities  of  transfer  as  a  function  of  v,  by  assuming 
rotational  equilibrium  in  DF(v)  and  CO^OO^O)  at  the  translation  temperature. 
We  assume  that  the  transfer  probability  decreases  about  one  order  of 
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magnitude  for  every  150  cm’'  of  energy  separation  as  Callear  (49)  found  for 
certain  diatomics.  In  the  final  states,  all  values  of  J  are  assumed  possible; 
however,  because  of  the  resonance  argument,  those  values  that  release  the 
smallest  quantity  of  energy  to  translation  are  most  likely.  At  300°K,  the 
higher  levels,  v  -  6,  7,  and  8,  are  still  most  likely  to  transfer.  A  calcula¬ 
tion  made  on  the  assumption  of  a  decrease  of  one  order  of  magnitude  for 
every  450  cm’'  of  energy  separation  gave  a  similar  result.  At  this  time,  a 
rate  of  transfer  based  on  these  resonance  and  equilibrium  population  argue- 
ments  seems  unwarranted,  especially  in  light  of  the  apparent  insensitivity  to 
resonance  exhibited  by  the  transfer  efficiencies  of  other  hydrogen  halides. 
Hence,  until  more  data  are  available,  we  estimate  the  rate  to  be  the  same 
for  all  v. 

We  use  the  values  of  Refs.  [35]  and  (47)  to  obtain 

k,0  =  4.  7  x  10^T"^ 

V 

Since  the  temperatures  used  to  determine  this  function  are  very  close 
together,  this  rate  coefficient  should  be  used  with  prudence  outside  the 
300°  to  350°K  temperature  range.  In  fact,  for  large  temperature  variations 
the  3.  82  x  10  cc/mole-sec  rate  of  Ref.  (35)  may  be  more  appropriate. 
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IV.  C02  KINETICS 


An  extensive  theoretical  study  of  the  rates  of  VV  and  VT  reactions  of 
C02-C02  and  C02-N2  has  been  made  by  Herzfeld  [50,  51],  We  use  his  work 
as  the  basis  for  selecting  the  reactions  shown  in  Fig.  1,  and  we  use  Ref.  [51] 
as  a  guide  for  estimating  rates  when  experimental  measurements  are  not 
available.  In  addition,  many  of  the  reactions  in  our  model  are  common  to 

1  l 

those  reviewed  by  Taylor  and  Bitterman  [5j.  The  states  (03  0)  and  (03  0)  are 
assumed  in  equilibrium  and  statistically  combined.  We  denote  this  combined 
state  as  (030)  and  adjust  the  kinetics  accordingly.  The  small  energy  difference. 
73  cm  *,  between  these  states  is  neglected,  and  both  states  are  assignee!  :he 
energy  of  (03*0). 

A.  VIBRATIONAL  -  VIBRATIONAL  TRANSFER  IN  C02 

Intermoleculaar  VV  transfer  processes  represent  important  and  fast 
mechanisms  for  energy  transfer  within  CO£.  These  reactions  have  a  signifi¬ 
cant  effect  on  the  relative  populations  of  the  C02  levels,  including  those  that 
participate  in  lasing. 

(11)  CO2(00°l)  +  C02(00°0)  CO2(10°0)  +  002(01*0)  +  294  cm”1 

(12f)  COJOO°1)  +  C0,(00°0)  ^CO,(02f0)  +  CO,(01*0)  cm-l  !  =  ? 

c  c  c  c.  cm  s  ~  c 

(13)  CO2(ll*0)  +  C02(00°0)  ^  CO2(10°0)  +  CO2(,01*0)  +  21  cm'* 
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(14)  C02(030)  +  C02(00°0)  5*  C02(10°0)  +  C02(01l0)  -  123  cm"1 

(15*)  C02(030)  +  C02(00°0)  ^  C02(02^0)  +  C02(01l0)  -{7° 

(16)  C02(10°0)  +  C02(00°0)  ^  CO2(01  X0)  +  CO2(01 10)  +  54  cm"1 

(17*)  C02(02*0)  +  C02(00°0)  ^  C02(0110)  +  C02(0110)  +  {“4J  ® 


Rate  coefficients  for  Reacts.  (11)  and  (12)  have  been  calculated  by 
Herzfeld  [5 1J  between  300°  and  3000° K.  We  have  approximated  these  rates 
as  a  function  of  T  to  be 


k„  .  34.  0  T2'  V-  3,0 


1  112.0  T2,  5e°  t  =  0 

K12  (neglected  1=2 

The  value  of  kj2  is  nearly  two  orders  of  magnitude  less  than  k^2  [51]. 

Theoretical  or  experimental  estimates  of  rates  of  Reacts.  (13)  -  (15)  are 

not  available.  There  is,  however,  a  laser  fluorescence  measurement  by  Rhodes, 

Kelly,  and  Javan  [52]  that  has  been  interpreted  as  either  the  rate  of  React.  (16)  or 

12 

(17).  At  296°K,  they  found  this  rate  coefficient  to  be  7.  5  x  10  cc/mole-sec. 

Sharma  [53],  using  the  theory  of  Sharma  and  Brau  [54,  55],  and  Seeber  [56],  using 

a  modified  SSH  theory,  have  calculated  the  rate  coefficient  of  React.  (16)  at 

12  12 

300° K  and  find  3.  2  x  10  and  2,5  x  10  cc/mole-sec,  respectively. 
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Sharma  [53]  suggests  that  the  experiment  of  Rhodes,  et  al.  [52]  is  controlled 
by  React.  (17  ),  which  he  finds  is  21  times  faster  than  React.  (16).  Seeber 
[56]  rejects  this  interpretation  and  suggests  that  the  measurement  reflects  a 
stimulated  Raman  effect,  multiple  photon  process,  and  a  subsequent  rotational 
relaxation.  Even  if  this  interpretation  is  correct,  this  effect  is  enhanced  by 
the  simultaneous  application  of  10.  6  and  9.  6-|i.m  radiation  in  the  experiment  of 
Ref.  [52]  and  would  not  arise  during  normal  lasing  conditions.  Nevertheless, 

interpretation  of  the  experiment  remains  in  question.  Seeber  [56]  calculates 

2  2  0 

%  8kj£,  and  Seeber  and  Sharma  both  find  k^  to  be  larger  them  k^  by  the 

factors  of  8  and  35,  respectively. 

Other  measurements  of  the  relaxation  of  the  (10^0)  state  [57  -  62J  give 
rate  coefficients  ranging  from  0,4  to  2.8  x  10^  cc/ mole-sec.  Interpretation 
of  these  measurements  is  very  difficult,  and  the  uncertainty  in  the  reported 
rate  coefficients  is  large. 

In  view  of  these  uncertainties,  we  estimate  the  rates  of  Reacts. (16)  and 
(17*)  from  the  SSH  calculations  of  Seeber  [56]  and  assume  that  the  probability 
is  independent  of  T;  hence,  we  set 


k16  =  1. 5  x  10nT0,  5 


1. 5  x  1011T0,  5 


Reaction  (13)  is  similar  to  React.  (17*)»  in  that  both  involve  the 
transfer  of  one  quantum  of  to  a  ground  state  molecule;  hence,  we  set 


1 


k 


13  = 


0.5  k 


0 

17 


Approximation  of  the  bending  mode  of  CO^  by  an  harmonic  oscillator  model 
suggests  that1 


0.75 


1.5 


Reaction  (14)  involves  the  transfer  of  one  v 2  quantum  and  the  excitation  of  Vj 

_2 

at  the  expense  of  Zv^.  We  estimate  that  this  process  is  10  times  as  efficient 
as  React.  (15),  i.  e.  ,  one  order  of  magnitude  slower  for  each  quantum  changing 
modes;  hence,  we  set 


k14  = 

B.  VIBRATIONAL  -  TRANSLATIONAL  TRANSFER 

We  discuss  two  types  of  reactions  in  this  section.  The  first  type  is 
the  VT  process  by  which  exactly  one  quantum  of  an  excited  mode  is 


l The  value  of  k^  is  used  as  a  reference,  since  the  degree  of  resonance  is 

2 

much  less  than  that  of  (17  )  for  these  reactions. 


-32- 


exchanged  for  translational  energy  during  a  collision.  An  example  of  such 
a  reaction  is 


CO2(0110)  +  M  ^  C02(00°0)  +  M  +  667  cm"1 

The  second  reaction  type  is  an  intramolecular  VV  process  that  leaves 
the  excited  molecule  in  a  different  mode  from  the  original  state.  These  pro¬ 
cesses  differ  from  the  intermolecular  VV  reactions  where  the  state  of  both 
molecules  is  changed  during  the  collision.  Again,  the  excess  energy  is 
exchanged  for  translational  energy.  An  example  of  such  a  reaction  is 

CO2(00°l)  +  M  ^  C02(1 1 1 0)  +  M  +  273  cm"1 

Both  processes  represent  transfer  of  vibrational  energy  to  translational 
energy  of  the  colliding  pair. 

Deactivation  of  the  upper  laser  level  by  collisions  represents  a  loss  of 
laser  energy,  whereas  deactivation  of  the  lower  laser  level  enhances  the 
inversion.  In  addition,  mechanisms  that  couple  levels  adjacent  to  and  below 
the  lasing  levels  are  important  for  realistic  computer  modeling.  In  Table  IV, 
Reacts. (18/  -  (29)  include  processes  of  this  type  considered  important  in  the 
modeling  of  C02  relaxation  for  transfer  laser  simulation.  We  first  discuss 
those  reactions  where  experimental  measurements  are  available,  and  then  we 
complete  the  discussion  by  making  theoretical  estimates  of  the  rates  of  the 
remaining  reactions. 
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1. 


DEACTIVATION  OF  THE  UPPER  LASER  LEVEL 


Reactions  (18)  and  (19)  represent  a  loss  of  laser  power  through 
collisional  deexcitation  of  00^(00^1) 


(18) 

CO2(00°l)  +  M 

^  CO2(ll10)  +  M  +  273  cm'1 

(19) 

CO2(00°l)  +  M 

^  C02(030)  +  M  +  417  cm"1 

Experimental  data  for  the  relaxation  rate  of  00^(00^1)  are  generally  pre¬ 
sented  without  specifying  the  final  state  of  CO^.  The  most  probable  final 

t 

states  are  those  of  Reacts.  (18)  and  (19).  One  might  also  expect  CO^(04  0) 
to  be  a  product  of  this  VT  reaction  of  00^(00^1 );  however,  we  have  limited 
our  consideration  to  levels  below  (00^1).  Herzfeld's  calculations  show  that 

3 

the  final  state  (03  0)  may  be  neglected,  while  the  ratio  for  pure  CO^  of  k^ 

(with  final  state  (03^0))  to  kjg  may  be  approximated  by  7.  7  x  10'^T®*  ®.  This 
ratio  is  0.  074  at  300°K  and  0.  19  at  1000°K;  hence,  React.  (18)  is  expected  to 
be  the  dominant  mechanism  over  this  temperature  range. 

Experimental  measurements  of  the  relaxation  of  mode  in  CO^  have 
been  made  by  several  techniques.  Laser  fluorescence  measurements  have 
been  made  by  several  investigators.  Hocker,  et  al.  [ 6 3 J  found  this  rate  to  be 
385  torr  * -sec  *  for  pure  C02  at  room  temperature.  Similar  measurements 
were  made  at  293®K  for  the  relaxation  of  by  collisions  with  He,  H^,  D^. 

H^O,  N^,  and  pure  C02  by  Moore,  et  al.[64]  and  by  Yardley  and  Moore  [65] 
for  collisions  with  He,  Ne,  Ar,  Kr,  Xe,  and  pure  C02>  Alexander,  et  al.  [66] 
measured  the  relaxation  rate  in  pure  C02  at  233°  and  293° K.  A  comprehensive 
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study  of  the  temperature  variation  of  this  relaxation  rate  in  pure  CO^  and 
has  been  made  by  Rosser,  Wood,  and  Gerry  [67J.  Rosser  and  Gerry  [68,  69 1 
have  also  measured  this  rate  over  a  wide  temperature  range  for  deactivation 
by  He,  O^,  H^O,  H NO,  and  Cl^.  The  temperature  dependence  of  relax¬ 
ation  for  deactivation  by  He,  Ne,  Ar,  Kr,  Xe,  HD,  and  over  the  range 
298°  to  800° K  is  given  by  Stephenson,  Wood,  and  Moore  [ 7 0 J.  Deactivation 
of  v ^  by  DF  and  HF  has  been  measured  by  Chang,  McFarlane,  and  Wolga  [ 7 1 J 
at  348°  and  373°  K  and  by  Stephens  and  Cool  [35J  at  350°  K.  With  the  use  of  a 
standard  modulated  light  source,  instead  of  a  laser,  Houghton  [72]  measured 
the  relaxation  of  in  pure  CO^  and  in  Ar.  Cheo  [61,  62]  measured  this 
relaxation  rate  at  300°Kfor  pure  CO^,  H^,  He,  N^,  and  other  species  with 
the  pulse-gain  method.  The  shock  tube  has  also  been  used  to  study  this 
relaxation  rate  by  Simpson  and  Chandler  [73]  for  mixtures  of  pure  CO^  and 
CC>2  with  Ar,  N^,  Ne,  He,  D^,  and  H^  over  a  wide  temperature  range.  In 
general,  their  results  are  in  poor  agreement  with  laser  fluorescence  measure¬ 
ments,  We  expect  that  the  process  involved  in  relaxation  of  in  the  shock 
tube  is  much  more  complex  than  that  in  the  laser  fluorescence  measurements, 
since  with  laser  excitation  only  one  state  (00^1)  is  produced  in  abundance  and 
interpretation  of  the  intramolecular  VV  data  is  more  straightforward.  The 
relaxation  times  observed  with  laser  fluorescence  are  essentially  a  direct 
measure  of  (00^1)  deactivation,  to  (11*0)  based  on  the  calculation  of  Herzfeld 
[51],  whereas,  in  the  shock  tube,  many  adjacent  levels  are  contributing  to  the 
activation  and  deactivation  of  ^  and  the  interpretation  of  pT  is  much  more 
complicated.  We  expect  that  detailed  modeling  of  the  CO^  kinetics  may  be 
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helpful  in  explaining  these  shock  tube  results.  Similar  problems  have  been 
found  for  measurements  made  with  the  spectrophone  technique  [74  -  79]. 
Levinson,  et  al.  [80]  have  measured  the  relaxation  rate  of  CO2(00°l)  by 
exciting  a  gas  mixture  with  a  current  pulse  and  monitoring  the  4.  3-^m  spon¬ 
taneous  decay.  Finally,  Biryukov,  Serikov,  and  Trekhov  [81]  measured  the 
rate  of  relaxation  of  induced  by  collisions  with  N^,  He,  and  H^O  and  in 
pure  CC>2  with  the  phase  method. 

Reduction  of  the  experimental  relaxation  times  to  rate  coefficients  is 
not  always  a  straightforward  process.  Taylor  and  Bitterman  [5]  have  used 
the  5 -reaction  model  of  Taylor,  et  al.  [82]  to  reduce  the  pT  data  for  (00  1). 
We  believe  that  it  is  just  as  appropriate  to  interpret  the  laser  fluorescence 
measurements  when  intermolecular  VV  processes  are  not  involved  as 
direct  measurements  of  kjg.  Hence,  we  use 


k 


18  = 


RT 

Ptco2-m 


(4) 


to  reduce  this  data.  The  ratio  of  k^  to  kjg  calculated  by  Herzfeld  is  then 
used  to  determine  kjg. 

Reduction  of  the  shock  tube  data  of  Ref.  [73]  represents  an  even  more 
difficult  problem.  We  have  not  used  these  data,  since  they  are  in  poor  agree¬ 
ment  with  other  measurements.  In  Fig.  5,  we  show  the  rate  coefficients  for 
the  various  species  of  interest  in  our  model.  For  modeling  simplicity,  rates 
for  all  species  have  been  approximated  by  one  of  three  different  functions. 
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For  CO^-CO^.  we  have  approximated  the  data  of  Rosier,  et  al.  [67]  and 
Rosser  and  Gerry  [69]  to  be 


k 


18a  " 


8.85  x  10‘4T4,75e“1,484° 


and  for  CO^-He,  we  find  that 


k18b 


=  1.13 


x  10-V*  8e"2*4360 


approximates  the  data  of  Rosser  and  Gerry  [68]  and  Stephenson,  et  al.  [70], 

1  5 

For  M  =  H2  (Fig.  5b),  we  find  that  a  T  *  temperature  dependence  matches 
the  data  of  Refs.  [69]  and  [70] .  Hence,  from  Refs.  [64]  and  [70],  we  estimate 
the  rate  coefficient  for  deactivation  by  to  be 

k18c  =  1.2*106t'-5 

In  addition,  we  estimate  from  Refs.  [35]  and  [71]  that  DF  is  4.5  times  as  effec¬ 
tive  as  H^  at  350° K  in  deactivating  vy  We  feel  that  the  large  CC^-DF  deac- 

1  5 

tivation  rate  at  350° K  indicates  a  weak  temperature  dependence  such  as  T  *  . 
Relaxation  of  by  Ar  is  approximately  as  fast  as  He  [70].  We  also  use  four 
times  the  CO^-He  rate  for  D  by  assuming  CC^-D^  and  CC^-D  are  about  the 
same  at  room  temperature,  and  we  estimate  the  rate  of  deactivation  by  F  to 
be  similar  to  DF.  We  assume  that  the  rate  coefficient  for  F^  is  similar  to  C^; 
hence,  we  take  CO^-F^  rate  to  be  twice  that  for  CO^-He.  In  Fig.  5a,  the  SSH 
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rate  for  C02_C02  calculated  by  Herzfeld  has  been  increased  by  a  factor  of 
10  for  graphic  comparison  with  experimental  data. 

2.  TRANSFER  BETWEEN  FERMI  RESONANCE  STATES 

Theoretical  SSH  calculations  by  Herzfeld  [51]  and  Seeber  [56]  for  pure 
CO2  indicate  that  the  rate  of  energy  transfer  between  the  strongly  coupled 
levels  in  Fermi  resonance 

(25f)  CO2(10°0)  +  M  ^  C02{02i0)  +  M  +  102.  8  cm-1  f=0 

is  fast.  However,  these  calculations  also  show  that  the  reaction  coupling  the 

2  0 

levels  not  in  Fermi  resonance,  (25  ),  is  faster  than  (25  ).  The  laser  fluor¬ 
escence  measurement  of  Rhodes,  Kelly,  and  Javan  [52]  at  300° K  in  pure  CO2 

indicates  that  k^  or  k^  is  possibly  10  times  faster  than  the  calculated  rate 

1  5 

coefficients  of  Ref.  [51],  which  vary  as  T  *  .  Since  interpretation  of  this 

experiment  remains  in  question,  we  use  the  results  of  Seeber  at  300°K,  and 

1  5 

the  rate  coefficients  are  assumed  proportional  to  T  *  , 

k®5  .  ax  io8t‘-5 

k*5  =  5.  8  x  108T1, 5 

Herzfeld's  calculations  also  show  that  ^  is  nearly  as  efficient  as  CO2  for 
these  reactions.  We  adjust  the  rate  coefficient  for  each  species  such  that  all 
have  nearly  the  same  probability  per  collision  as  CO Note  the  rate  of 
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Reacts,  (25*)  during  lasing  will  always  be  less  than  the  transfer  rate  of  DF 
to  CO^. 

The  transfer  between  the  other  states  in  Fermi  resonance,  (11*0,  03*0), 
is  given  in  React.  (20): 


(20)  ^(ll1©)  +  M  C02(03*0)  +  M  +  144  cm'* 

Although  Herzfeld  does  not  compute  this  rate,  he  does  compute  the  total  rate 
of  deactivation  of  the  combined  (11*0,  03*0)  state  to  be  20  times  slower  than 
React.  (25^)  at  300°  K.  This  is  probably  a  lower  limit  on  Ic^q*  we  e8t*mate 


no 


In  our  model,  React.  (20)  includes  the  relaxation  between  the  levels  (11*0) 

3 

and  (03  0)  that  are  not  in  Fermi  resonance.  In  any  case,  React.  (20)  is  of 
very  little  importance  during  lasing. 

3.  DEACTIVATION  OF  THE  BENDING  MODE  OF  COz 

Relaxation  of  mode  v ^  of  CO2  has  been  widely  studied  (83  -  86).  Taylor 
and  Bitterman  1 5)  have  compiled  relaxation  rates  fr  r  this  reaction 

(29)  CO2(01  *0)  +  M  ^  C02(00°0)  +  M  +  667  cm'* 


and  they  indicate  a  "best  fit"  over  the  temperature  range  of  interest  for 
M  =  CO2,  He,  H2  and  for  other  species  not  of  interest  in  our  model. 
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Conversion  of  the  recommendations  of  Ref.  [  5)  to  forward  rate 


coefficients  with  Eq.  (1)  yields  for  M  =  CO^ 


k29a  =  2‘  17  *  10"2t4'  19e’1#  1320 


For  He,  we  use 


k29b  =  8  27  X  102T 


2.95  -0.84336 
e 


The  recommended  value  for  M  =  decreases  with  T.  Sharma  [85],  using 
the  theory  of  Sharma  and  Brau  (54,  55]  ,  has  found  the  efficiency  of  in 
deactivating  v ^  to  vary  approximately  as  T~^,  and  the  agreement  of  theory 
and  experiment  is  good.  Thus,  for  M  =  D^,  we  assume  the  rate  is  propor¬ 
tional  to  T"1  and  adjust  it  to  the  value  reported  by  Cottrell  and  Day  [86]  at 
303° K,  since  their  values  for  M  =  C02  and  He  are  close  to  the  "best  fit" 
curves  of  Ref.  [5].  In  addition,  Simpson  and  Chandler  [73]  have  determined 
this  rate  from  their  shock  tube  measurements.  Although  these  values  are  in 
closer  agreement  with  other  measurements  than  their  observed  relaxation 
curves  for  (OO01),  we  have  not  used  their  values  in  the  present  model.  The 
temperature  variation  of  these  rate  coefficients  is  shown  in  Fig.  6,  along 
with  values  calculated  by  Herzfeld  (51)  for  pure  C02.  The  overall  agreement 
of  the  theoretical  and  experimental  rates  is  much  better  than  that  found  for 
relaxation  of  vy 

Relaxation  times  for  C02>Ar  have  been  reported  by  Huety  and 
Chevalier  [84]  and  Simpson  and  Chandler  [73].  The  impact  tube  measurements 
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RATE  COEFFICIENT  (cc/mole -sec) 


T  (*K) 


Fig.  6.  Experimentally  derived  rate  coefficients  for 
deactivation  of  COgfOl^O)  by  M:  M  =  COg 
and  He,  Ref.  [5j;  O  Ref.  [84],  Q  Ref.  [86]; 
X  SSH  calculations,  Ref.  [51 J. 
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of  Ref.  I  84)  are  more  consistent  with  data  of  other  investigators;  thus,  we 
take  their  results  and  assume  Ar  to  be  6%  as  efficient  as  CO^  in  deactivating 
v<2*  Relaxation  of  i- 2  by  F,  F^,  DF,  and  D  has  not  been  measured.  The 
efficiencies  of  these  species  are  estimated  to  be  similar  to  C02  with  the 
exception  of  D,  which  is  assumed  similar  to  D^.  These  arc  gross  approxima 
tions,  and  additional  measurements  for  these  species  are  needed. 

4.  OTHER  REACTIONS 

The  rates  of  the  remaining  reactions  have  not  been  measured.  These 
reactions  include 


(21*) 

CO2(ll10)  +  M 

^  CO,(02f0)  +  M  +  [ZJj  "m- 

2'  |741  cm 

(22) 

CO^l^O)  +  M 

^  CO2(10°0)  +  M  +  688  cm"1 

(23*) 

C02(030)  +  M 

^  CO, (02*0)  +  M  +jj?^  Cm"1 

2  (597  cm 

(24) 

C02(030)  +  M 

^  CO2(10°0)  +  M  +  544  cm"1 

(26) 

CO2(10°0)  +  M 

^  CO2(0110)  +  M  +  721  cm"1 

(27) 

C02(0220)  +  M 

^  C02(02°0)  +  M  +  50  cm"1 

(28*) 

C02(02*0)  +  M 

-co2(oi‘o)  +  M  +  j^  l<: 
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Theoretical  estimates  of  the  rates  of  these  reactions  have  been  made  by 
Herzfeld[5l]  for  M  =  CO,.  Our  analytical  approximations  of  his  r  \tes  are 


k®,  =  S.ixlO-V-'V-9030 


k*, 


k„  =  5.  9T3'  78e°-  5490 


k°3  =  3. 87  x  102T3,3e1,24° 


k23  =  5.7  x  102T3,  32e1,22° 


k24  =4.9xl02T3*04e1*060 


k26a  =  4‘  4  v  102T3,26e1,4f 


K„‘  5.9x108T1-5 


.0  _  ,  ,  v  lft3T3.  1  1.40 

k28a  “  2*4xl°  T  e 


.2  _  ,  1n2_3.  31  l-  30 

k28a  =  6.0x  10  T  e 


We  assume  both  levels  of  vibrational  angular  momentum  in  (030)  are 
deactivated  at  the  same  rate  in  determining  k^3  and  k24< 


Seeber  [56]  has  also  calculated  at  300*  K,  and  his  value  is  35%  larger 
than  that  of  Herzfeld  and  has  been  used  with  the  temperature  variation  calcu¬ 
lated  by  Ref.  ( 5 1 1  to  determine  k^y 

Reactions  (21)  through  (24)  represent  collisional  mechanisms  that  may  con¬ 
tribute  to  the  population  of  the  lower  laser  level  from  products  of  deexcitation 
of  the  upper  laser  level.  Since  experimental  estimates  of  the  rate  of  these 
reactions  are  not  available,  we  use  Herzfeld's  rate  coefficients  and  assume 
the  probability  for  reaction  is  the  same  for  each  species.  In  addition,  we 
assume  that  all  species  are  equally  efficient  in  equilibration  of  the  vibra¬ 
tional  momentum  of  (02*0),  React.  (27). 

For  Reacts.  (26)  and  (28),  we  use  the  theoretical  rate  coefficients  of 
Herzfeld  for  species  (Table  IV).  Since  the  temperature  variation  of 

React.  (29)  for  an<*  M  ^  is  considerably  different  from  the  calculation 
of  Herzfeld,  we  use  the  temperature  dependence  of  the  rates  of  React.  (29) 
for  these  species  and  adjust  the  relative  rates  of  the  reaction  to  be  in  the 
same  ratio  as  that  calculated  by  Herzfeld  for  CO^-CO^  at  1000°K,  where  the 
relative  magnitude  of  the  calculated  rates  may  be  expected  to  be  fairly  accu¬ 
rate,  At  1000#K  for  CO-j-CO^,  Herzfeld  finds 

k26:k28:k28:k29  =  8:3.  4:4*  0:  0 
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Hence,  we  take 


k26b  =  I,8k29b 
k26c  =  1,8k29c 

k28b  =  3,4k29b 


and  continue. 

The  complete  kinetic  model  for  the  DF-CO^  chemical  transfer  laser  is 
given  in  Table  IV. 


V.  COMPUTER  EXPERIMENTS 


A.  LASER  SIMULATION 

The  chemical  laser  computer  simulation  used  in  this  study  is  described 
in  Ref.  [2] .  A  resume  of  its  features  is  presented  here.  Rate  equations  are 
used  to  represent  the  chemical  kinetic  and  stimulated  emission  processes 
occurring  in  a  representative  unit  volume  within  a  Fabry-Perot  cavity.  All 
processes  are  assumed  to  be  uniform  throughout  the  cavity.  Only  the  (00%)  - 
(10%)  band  is  permitted  to  lase,  and  a  Boltzmann  distribution  of  the  rotational 
levels  at  the  translational  temperature  T  is  assumed.  Only  the  transition 
with  maximum  gain  is  assumed  to  lase,  and  this  is  always  in  the  P-branch. 
During  lasing,  the  gain  at  line  center  is  held  constant  at  its  threshold  value 

O..  • 

thr 

The  chemical  reactions  may  be  written  as 


Z>riN  k*I>riN  l5» 

i  -r  i 

where  |Nj|  is  the  molar  concentration  of  species  i,  and  are  stoichio¬ 
metric  coefficients,  and  k±r  are  the  forward  and  backward  rate  coefficients. 
The  rate  of  change  of  concentration  of  CO2(00%)  is  given  by 


d[cO2(00°l)| 


dt 


-Xrad<J>  +  *ch<00 


(6a) 
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For  CO2(10 


u),  we  have 


d|c02(10°0)| 

dt 


+  Xch'l()0<» 


(6b) 


and  for  all  other  species 

d  In.  I 

dT=xch<i)  (6‘> 

The  photon  emission  rate  Xracj(J)  is  the  rate  of  change  of  concentration  result¬ 
ing  from  lasing,  and  J  is  the  rotational  quantum  number  for  the  lower  level 
of  the  maximum  gain  transition  within  the  10.6^im  band.  The  rate  of  change  of 
species  concentration  resulting  from  chemical  reactions  is 


where 


xch<‘>  ■  I>ri  -  %i>Lr 


(7) 


The  laser  cavity  is  assumed  to  have  a  uniform  photon  flux  with  active 
medium  length  L  and  mirror  reflectivities  Rq  and  R^.  Lasing  initiates  when 
the  highest  gain  transition  of  the  .band  is  equal  to  the  threshold  gain  a^r  where 

*thr  ■  -2Tln<R0RL>  (8> 
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Only  P-branch  transitions  need  be  considered;  thus,  the  gain  of  a  transition 
with  lower  level  J  is 


«(J)  3  -^rw(J)B(J)<MJ)||f^-}lCO2(000l,J-l)]  -  lCO2(10°0,J)jj  (9) 

where  the  wave  number  of  the  transition  is  w(J)  and  B(J)  is  the  Einstein 
isotropic  absorption  coefficient  based  on  the  intensity.  Line-broadening  con¬ 
stants  and  resonance  constants  used  in  the  Voigt  profile  $(J)  at  line  center 
are  those  of  Ref.  [2] . 

The  Boltzmann  distribution  of  the  rotational  populations  is  given  by 

[C02(v1(w£,  v3;J)]  =  [C02(^,wj,  w3)]  e‘hcEJ/kT  (10) 

r 

where  t/20^  is  the  rotational  partition  function,  the  additional  factor  of  2  is 
a  statistical  symmetry  factor,  Gy  is  the  characteristic  rotational  tempera¬ 
ture  [2],  and  rotational  energy  Ej  is  from  Herzberg  [20] .  We  note  that  only 
odd  J  values  are  present  in  the  upper  states  and  even  values  in  the  lower 
states.  Planck's  constant,  the  speed  of  light,  and  Boltzmann's  constant  are 
denoted  as  h,  c,  and  k,  respectively.  The  energy  equation  for  a  constant 
density  gas  is  written  in  the  form 
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where  is  the  specific  heat  at  constant  pressure.  H.  is  the  molar 
enthalpy  of  species  i.  p  is  the  pressure,  and  is  the  output  lasing  power 
per  unit  volume  given  by 

PL(t)  *  hcNAu>(J)Xrad(J)  (12) 

The  average  intensity  inside  the  cavity  may  be  calculated  from  the 
relation 

I  =  PL(t)L/(l  - 

A  fourth-order  Runge-Kutta  routine  is  used  to  solve  Eqs.  (6),  (11),  and 
the  equation  of  state  to  determine  pressure,  temperature,  and  species  con¬ 
centrations  until,  at  some  time  tg  and  for  some  value  J  =  Jq,  the  gain  on  a 
given  vibration- rotation  transition  reaches  At  this  time,  the  laser  pulse 

begins.  Then  Eqs.  (6),  (8),  (11),  (12),  and  the  equation  of  state  are  solved 
for  the  transient  temperature,  pressure,  concentrations,  output  power,  and 
active  J.  During  lasing,  this  value  of  J  shifts  as  a  result  of  the  temperature 
rise  and  concentration  changes.  Lasing  terminates  when  all  gains  drop  below 

°thr  ’ 

The  pulse  energy  per  unit  volume  is  given  by 

E  =  /t0+tcPLdt  (13) 

*0 
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where  t  is  the  pulse  duration.  For  simplicity,  initiation  is  modeled  by  the 
introduction  of  a  finite  concentration  of  F  atoms  into  the  initial  gas  mixture. 
B.  CHEMICAL.  COLLISIONAL.  AND  CAVITY  MECHANISMS 

The  relative  importances  of  the  various  kinetic  processes  are  examined 
by  simulating  laser  performance  at  typical  laboratory  conditions.  Two  initial 
gas  mixtures  of  XF:lF2:lD2:40He:8CO2  at  300*K  and  50  torr  with  X  3  0. 1 
and  0.01  are  used,  where  the  F  atoms  represent  varying  degrees  of  dis- 
sociation  of  by  some  initiation  method.  Mirror  reflectivities  are  0.9  and 
1.0,  and  the  active  medium  length  is  53.3  cm.  The  following  paragraphs  give 
an  interpretation  of  the  principal  kinetic  mechanisms  of  these  cases. 

By  noting  the  conservation  of  vibrational  quanta  in  DF(v)  and  properly 
accounting  for  the  vibrational  levels,  one  may  deduce  the  relation 

£VPv  =£Tv*E,™*^Dv1*-I  (14) 

v*l  v*l  v»l  v«l 

where  is  the  rate  of  pumping  into  level  v  by  Reacts.  (4)  and  (5),  Ty  is 
the  rate  of  transfer  from  DF(v)  to  CO2(00®l)  by  React.  (10),  dDF(v)/dt  is  the 
rate  of  change  of  the  DF(v)  concentration,  and  D  ,  is  the  total  rate  of 

V  f  Vm  X 
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deactivation  of  DF(v)  to  DF(v-l)  by  Reacts.  (6)  and  (8).^  The  terms  in 
Eq.  (14)  may  be  interpreted  as 


Chemical 

Collisional 

pumping 

= 

transfer 

rate 

rate 

Level  filling 
rate  in  DF 


+ 


Collisional 
deactivation 
rate  of  DF 


The  rate  of  excitation  of  CC^OO^l)  by  DF(v)  is  equal  to  the  transfer  rate; 
therefore,  conservation  of  number  density  of  CC^OO^l)  may  be  expressed  as 


(15) 


where  is  the  rate  of  collisional  deactivation  of  (00  1). 

The  relative  importance  of  these  mechanisms  throughout  the  pulse  is 
graphically  illustrated  in  Fig.  7.  An  abbreviated  notation  for  the  terms  in 
Eqs.  (14)  and  (15)  is  used  where  the  contribution  from  the  pumping  React.  (4) 
for  v  =  1,  •  •  *  ,  4  is  referred  to  as  P^,.  For  v  >  5,  React.  (4)  normally 
transfers  matter  in  the  reverse  direction  and  is  denoted  as  P^.  The  contribu¬ 
tion  from  React.  (5)  is  labelled  P^.  The  total  rate  of  collisional  deactivation 
of  DF(v)  is  denoted  by  D^.  The  total  collisional  transfer  rate  is  T,  and  the 
rates  of  level  filling  for  DF(v)  and  COgtOO^l)  are  denoted  as  Lp  and  L^., 
respectively. 


^The  rate  of  React*  (3)  is  normally  too  small  to  contribute  to  this  term,  and 
the  number  of  quanta  is  not  changed  by  the  VV  processes,  React.  (7). 
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IMTIU.  CONDITIONS 

C4NITY 

F  f,  Dj  COj  Ht 

l  >  553t« 

001  1  1  1  40 

V  10 

P,  50  Tor 

»0‘  0’ 

T,<  300'« 

o  So  tt>  iso  wo  no 

TIME  (pwc) 


Fig.  7.  Relationship  of  dominant  kinetic 
mechanisms  in  the  DF-CO2 
model,  (a)  Initial  F  =  0.  1. 

(b)  Initial  F  =  0.01. 


We  examine  the  mechanisms  for  F  3  0.1,  Fig.  7a,  in  detail  and  point 
out  only  those  features  that  are  different  for  F  3  0.01  (Fig.  7b).  Approxi¬ 
mately  60%  of  the  initial  F  atom  concentration  is  consumed  in  the  first  few 
microseconds  by  React.  (4),  P  .  During  this  period,  D  atoms  are  produced, 
and  the  vibrational  levels  of  DF(v)  begin  filling.  During  this  transient,  the 
rate  of  React.  (5) increases  with  onset  of  operation  of  the  chain,  and  transfer  to 
CO^tOO^l)  is  followed  by  emission  of  lasing  photons  when  cavity  threshold  is 
reached.  After  the  large  initial  consumption  of  F  atoms,  the  concentration  of 
F  atoms  increases  slowly  and  monotonically,  and  we  have  quasisteady  chain 
operation  such  that 

1*1, =0  ■  1*1, ,o  ♦  M.,0 

Time  histories  of  species  concentrations  for  initial  F  =  0. 1  are  shown  in 
Fig.  8. 

The  rate  jf  level  filling  in  DF,  L^,  represents  a  rate  of  energy  storage. 
Some  of  this  stored  energy  is  lost  to  collisional  deactivation  at  a  later  time. 
From  Fig.  7a,  we  see  that  only  a  small  quantity  of  the  stored  energy  in  DF 
is  available  at  pulse  termination.  A  small  rate  of  level  filling  in  CO^tOO^l), 
L^,  is  required  to  maintain  the  necessary  gain,  Eq.  (8).  This  rate  is  cou¬ 
pled  to  the  rate  of  deactivation  of  the  lower  level,  Xra^(J)»  D^,  and  the 
transfer  rate. 

Lasing  continues  until  the  rate  of  transfer  into  ( 00®  1 )  is  not  sufficient 
to  offset  the  loss  mechanisms,  and  all  o(J)  drop  below  Pulse  termina¬ 

tion  is  caused  in  part  by  three  processes.  First,  depletion  of  the  fuel 
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supply  (D^  and  F^)  results  in  a  lower  mass  transfer  rate  for  Reacts.  (4)  and 
(5);  this  is  not  the  case  for  F  *  0.01.  Second,  deactivation  of  DF(v)  and  of 
CO^OO^l)  by  Reacts.  (18)  and  (19)  increases  sharply  near  pulse  termination 
because  of  the  increasing  DF  and  00^(00^1)  concentrations.  Recall  that  DF 
is  a  very  efficient  collisional  deactivator  for  both  DF(v)  and  00^(00^1).  Third, 
the  temperature  rise  resulting  from  the  exothermic  chain  reaction  gives  rise 
to  higher  concentrations  in  the  (01*0)  and  (02*0)  levels  of  CO^;  therefore, 
depletion  of  the  lower  laser  level  is  less  efficient. 

The  principal  differences  in  mechanisms  for  the  two  values  of  initiation 
are  easily  seen  in  Figs.  7a  and  7b.  With  the  higher  pumping  rate  for  F  3  0.1, 
the  transfer  rate  is  not  fast  enough  to  deplete  the  excited  levels  of  DF;  hence, 
the  rate  of  level  filling  is  higher  and  a  larger  DF(v)  population  results  in  higher 
deactivation  losses  in  DF.  Pumping  for  F  3  0.01  is  slower  and  is  still 
increasing  when  the  pulse  terminates.  Although  transfer  to  CO-,  is  more 
efficient  with  this  lower  pumping  rate,  only  50%  of  the  fuel  is  consumed  before 
pulse  termination,  as  compared  to  80%  for  F  3  0.1.  Although  deactivation  in 
CO^  is  more  important  for  F  =  0.  01  and  is  the  cause  of  pulse  termi- 

3 

nation,  this  is  partly  offset  by  the  lower  rate  of  deactivation  of  DF.  Hence, 
even  with  much  less  of  the  fuel  consumed  during  the  pulse,  the  chemical 
efficiency  for  F  =  0.  01  is  6%  as  compared  to  8%  for  F  =  0.  01.  A  more 
detailed  investigation  of  these  mechanisms  is  now  given. 


Lower  levels  of  initiation  are  often  found  to  give  rise  to  increased  importance 
of  mechanisms  deactivating  the  lasing  species  143]. 
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in  Fig.  9,  the  time  histories  of  the  rate  Lr  of  DF(v)  deactivation 
reactions  are  shown.  For  F  =  0.  1,  the  rapid  production  of  DF(v)  in  the 
upper  vibrational  levels  after  initiation  leads  to  deactivation  of  DF  by  the 
DF-D^  VV  React.  (8).  This  rate  eventually  decreases  as  the  concentration 
of  D^(v)  (v  >0)  increases  relative  to  DF(v),  until  the  reaction  begins  to  run 
in  reverse  near  termination,  as  shown  by  the  dashed  curve.  The  principal 
catalysts  deactivating  DF  are  F  and  D,  React.  (6e),  and  self- relaxation 
of  DF,  React.  (6c).  Catalysts  of  other  VT  mechanisms  included  in  the 
figure  are  CC^,  React.  (6a),  He,  React.  (6b);  and  D React.  (6d). 

A  lower  overall  rate  of  DF  deactivation  is  found  when  initial  F  =  0.01, 
and  React.  (6e)  is  less  significant  as  a  result  of  the  lower  F  and  D  concen¬ 
trations.  A  slower  DF-D^  VV  transfer  rate,  React.  (8),  is  also  found. 

With  initial  F  =  0. 1,  the  dominant  deactivation  catalyst  for  CC^OO0!) 
is  DF,  React.  (18c),  in  Fig.  10.  The  catalyst  for  this  reaction  also  contains 
F  and  D^i  however,  the  contribution  of  these  species  is  at  least  one  order 
of  magnitude  less  than  that  of  DF  for  most  of  the  pulse  duration.  The  high 
concentrations  of  He  and  CC^  cause  significant  deactivation  by  Reacts.  (18b) 
and  (18a).  As  a  result  of  our  earlier  assumption,  the  relative  rates  of 
React.  (19)  are  very  similar  to  those  of  (18).  These  deactivation  mechanisms 
are  very  similar  for  F  =  0.01,  although  the  rate  is  lower  as  a  result  of  the 
slower  transfer  rate  and  a  lower  concentration  in  (00^1). 

The  principal  mechanism  populating  the  lower  lasing  level  (10®0)  is 
laser  action,  *ra(j(J)»  Fig.  11.  Reactions  (1 3),  (22),  and  (24)  deactivate 
(11*0)  and  (030);  these  levels  are  supplied  by  deactivation  of  (00^1)  through 
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COLLI SIONAL  DEACTIVATION  RATE  Lr  (mole/cc- 


TIME  (psec) 

Fig.  10.  Time  history  of  the  rate  Lr  of  reactions  deactivating  the 
upper  laser  level. 
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.  Kinetic  mechanisms  that  produce  CC^UO^O). 
The  time  of  a  shift  in  the  lasing  transition  to 
a  new  value  of  J  is  labeled  with  the  new  J. 

(a)  Initial  F  =  0.1.  (b)  Initial  F  =  0.01. 
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Reacts.  (18)  and  (19).  The  increased  rate  of  React.  (13)  relative  to 
Reacts.  (22)  and  (24)  for  F  =  0. 01 ,  as  compared  to  F  a  0. 1 ,  is  a  result 
of  the  slower  rate  of  transfer  to  (00^1)  and  thus  a  lower  concentration  in 
this  level  and  a  higher  concentration  in  the  ground  state. 

Collisional  processes  that  deactivate  the  lower  laser  level  for  F  =  0. 1, 
Fig.  12,  and  for  F  s  0,01  are  similar.  Deactivation  of  ( 10^0)  to  (02*0) 
through  React.  (25*)  is  the  major  mechanism  that  empties  the  lower  level. 
The  intramolecular  VV  reaction  (1 0°0)-*(01  * 0),  React.  (26),  is  also  an 
important  process,  especially  for  M  =  He,  React.  (26b).  The  decrease 
during  the  reaction  in  the  efficiency  of  the  intermolecular  VV  process, 

React.  (16),  is  a  result  of  the  increase  in  temperature  and  equilibrium  popu- 
lation  of  (01  *0). 

Since  deactivation  of  (10°0)  is  closely  linked  to  the  (02*0)  levels  which 
in  turn  deactivate  through  (01*0)  to  the  ground  state,  we  examine  those  pro¬ 
cesses  that  depopulate  (02*0)  in  Fig.  13  for  F  =  0.1.  Again,  He  is  the 
principal  deactivation  catalyst,  React.  (28  *b).  Other  intramolecular  VV 
catalysts  CC^,  Reacts  (28 *a),  and  and  D,  React.  (28 *c),  contribute  very 
little  to  the  equilibration  of  these  levels.  Reactions  (17*)  are  significant 
early  in  the  pulse  as  VV  deactivation  processes,  but  near  termination  the 
increased  population  of  (01*0)  and  the  higher  temperature  cause  these  reac¬ 
tions  to  transfer  matter  in  the  reverse  direction.  The  relative  rates  of  these 
reactions  for  initial  F  -  0.01  are  similar  to  those  shown  in  Fig.  13. 

In  the  present  model,  deactivation  rate  coefficients  for  (02*0)  by  VT 
and  intramolecular  VV  processes  are  based  on  measurements  of  VT 
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Fig.  13.  Collisional  processes  that  equilibriate  (02*0). 


deactivation  of  (01*0).  Thereforei  the  VT  mechanisms  that  equilibrate 
(01*0),  React.  (29),  are  comparable  in  importance  to  the  corresponding  VT 
mechanisms  that  equilibrate  (02  *0),  React.  (28*).  For  the  mixtures  con¬ 
sidered  in  this  study,  we  also  find  that  state  (10^0)  and  all  levels  below  it  are 
very  near  equilibrium  throughout  the  reaction. 

C.  EFFECT  OF  VARIATIONS  IN  UNCERTAIN  RATE  COEFFICIENTS 

The  variations  in  rate  coefficients  from  the  standard  are  selected  to 
include  the  major  uncertainties,  at  the  time  of  this  study  (early  1972),  in 
rates  of  dominant  processes.  Of  course,  the  degree  of  uncertainty  of  a 
rate  coefficient  is  only  speculative  and  significant  changes  in  the  initial  mix¬ 
ture  composition  can  change  the  relative  importance  of  dominant  mechanisms. 
Hence,  the  following  variations  only  give  an  indication  of  the  sensitivity  of 
laser  performance  to  some  of  the  more  important  mechanisms.  The  effect 

of  these  changes  on  the  peak  power  Pl_.vi  pulse  energy,  and  pulse  duration 

max 

are  shown  in  Table  V. 

The  variations  from  the  standard  model  (Table  IV)  studied  and  their 
principal  effects  are: 

Rate  1  The  rate  of  DF-CO^  transfer,  React.  (10),  is  assumed  con- 

1 2 

stant  in  temperature  and  equal  to  3.  82  ^  10  cc/mole-sec  for 
all  DF(v).  There  is  some  evidence  that  the  negative  tempera¬ 
ture  dependence  in  the  standard  model  is  too  strong  [87].  When 
the  pumping  rate  is  large  relative  to  the  transfer  rate,  larger 
deactivation  losses  occur  before  transfer  from  DF(v)  to  CO^ 
takes  place.  With  kjQ  assumed  constant,  deactivation  losses 
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Table  V.  Effect  of  Rate  Coefficient  Changes  on  Pulse  Characteristics 


in  DF  are  reduced.  Hence,  this  change  gives  higher  values 

for  the  pulse  energy  and  peak  power,  and  the  effect  is  larger 

for  F  =  0.  1  where  the  pumping  rate  is  higher. 

Rate  2  Transfer  for  DF(1)  is  the  same  as  for  Rate  1  and  kjQv  =  vkiOj 

Preliminary  experimental  evidence  indicates  that  the  rate  of 

transfer  may  increase  with  v  [87]  .  By  comparing  these 

results  with  those  of  Rate  1,  both  Pr  and  E  increase  with 

"max 

this  change,  especially  for  F  =  0.  1.  The  competition  between 
deactivation  in  DF  and  transfer  to  CO^  indicated  earlier  is 

responsible  for  this  behavior. 

Rate  3  The  efficiency  of  deactivation  of  (00®  1)  by  and  F,  React.  (18) 
and  (19).  is  assumed  the  same  as  CO^-CO^.  The  laser  pulse 
characteristics  are  insensitive  to  thes<j  reactions  for  our  choice 
of  initial  concentrations. 

Rate  4  Rate  3  is  modified  such  that  the  efficiency  of  DF  in  deactivating 

(00®  1)  is  proportional  to  T"^  as  indicated  by  combining  results 

of  Ref.  [35]  with  preliminary  results  of  Cool  at  300°  K  [47],  Of 

course,  this  change  gives  rise  to  an  increase  in  E  and  Pl_ 

max 

since  deactivation  is  important  later  in  the  reaction  when  the 
temperature  and  the  DF  concentration  are  both  higher. 

Rate  5  The  rate  coefficients  for  all  processes  that  deactivate  (10®0) 
are  arbitrarily  reduced  by  a  factor  of  10.  This  decrease  was 
not  sufficient  to  cause  (10®0)  to  deviate  from  near  equilibrium; 
therefore,  the  pulse  characteristics  are  essentially  the  same 
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Rate  6 


Rate  7 


Rate  8 


Rate  9 


as  those  computed  with  the  standard  rate.  A  much  larger 
decrease  in  these  rate  coefficients  is  necessary  before  pulse 
characteristics  are  substantially  altered. 

Rate  coefficients  for  Reacts.  (13),  (22),  and  (24)  are  increased 
by  a  factor  of  10.  This  change  is  not  large  enough  to  prevent 
the  reactions  that  relax  (10°0)  from  maintaining  the  near 
equilibrium  condition;  thus,  no  change  is  found  in  pulse 
characteristics. 

Rate  coefficients  for  all  processes  in  Reacts.  (18)  and  (19)  are 

arbitrarily  increased  by  a  factor  of  2.  This  pessimistic  change 

is  greater  than  the  uncertainty  in  the  rates  of  these  processes 

and  causes  a  reduction  in  Pi  ,  E,  and  t_  +  t  . 

■“max  U  c 

The  overall  pumping  rates  and  k,.  are  increased  by  a  factor 
of  2.  As  expected,  ^Lmax  *ncreases  an<*  the  pulse  duration 
decreases  with  the  faster  chain.  Again,  the  competition 
between  deactivation  in  DF  and  transfer  to  CO^  manifests 
itself,  and  no  increase  in  E  for  F  -  0.  1  is  shown  while  E 
increases  for  F  =  0.  01  where  this  competition  is  less. 

The  distribution  of  D  +  F^  pumping.  Reacts.  (5),  are  assumed 
to  be  more  sharply  peaked  near  v  =  8.  The  new  distribution 
is  obtained  from  Herbelin's  new  interpretation  [88]  of  recent 
H  +  F^  data  of  Jonathan  and  co-workers  on  the  basis  of  revised 
estimates  of  the  HF  matrix  elements.  The  altered  distribution  is 
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Rate  1 


Rate  11 


Rate  12 


Rate  13 


k5(0):  k5(l):  kg{2):  k5<3):  k5(4):  k5<5):  k&(6):  k5(7):k5<8):k5(9)  = 
0:0.  052:0.  065:0.  078:0.  078:0.  208:0.  494:0.  792:  1.  0:0.  701 
With  this  distribution.  React.  (5)  converts  53%  of  the  available 
chemical  energy  into  vibrational  energy;  the  standard  rate 
converts  46%  of  this  energy.  The  values  of  Pi  ,  E,  and 
tg  +  tc  all  increase  with  this  change  as  expected. 

All  DF-DF  VV  rate  coefficients  are  taken  as  l/4th  the 
binary  collision  rate.  Preliminary  experimental  measurements 
by  Bott  (87]  indicate  that  these  rate  coefficients  may  be  this 
large.  As  in  the  HF  system  [43,44].  laser  performance  is 
insensitive  to  these  processes. 

The  rate  coefficient  for  deactivation  of  DF(v)  by  D  and  F  is 
increased  by  a  factor  of  10.  The  effect  of  this  rate  increase 
is  more  important  for  F  =  0. 1  where  the  F  and  D  concentra¬ 
tions  are  larger. 

The  efficiency  of  CO2  in  deactivating  DF  is  increased  by  a 
factor  of  10.  For  these  mixture  compositions,  this  change 
produces  a  very  small  decrease  in  the  laser  performance 
characteristics. 

The  DF-D^  VV  rate  coefficients,  React.  (8),  are  increased 
by  a  factor  of  10.  Laser  performance  for  these  concentra¬ 
tions  is  insensitive  to  these  reactions. 
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VI.  SUMMARY  AND  CONCLUSIONS 


A  kinetic  model  of  a  DF-CO^  chemical  transfer  laser  and  a  review  of 
the  rate  coefficient  literature  needed  to  formulate  this  model  are  presented. 

The  portion  of  this  model  is  applicable  for  modeling  the  DF  chemical 

laser.  Rate  coefficients  recommended  for  the  model  are  listed  in  Table  IV. 

A  computer  simulation  is  used  to  examine  the  relative  importance  of  various 
kinetic  mechanisms  during  lasing.  Two  initial  mixture  compositions  of 
XF  :  1  F^  :  1  :  8CO2  :  40  He  at  50  torr  and  300°  K,  with  initiation  levels 

modeled  by  setting  X  =  0.  1  and  0.  01,  are  studied.  For  F  =  0.  01,  a  260-psec 
pulse  with  6%  chemical  efficiency  is  predicted,  while  the  higher  level  of  initia¬ 
tion  produces  a  52-psec  pulse  with  8%  chemical  efficiency.  To  the  best  of  our 
knowledge,  this  is  the  first  published  calculation  of  a  CO^  lasing  system  that 
includes  detailed  modeling  of  the  rate  processes  that  relax  the  lower  laser 
level.  For  the  initial  compositions  considered,  the  (10^0)  level  and  those  below 
it  are  near  equilibrium.  Calculations  with  other  compositions,  however,  indi¬ 
cate  that  this  is  not  always  true  and  detailed  modeling  is  sometimes  necessary. 

The  sensitivity  of  laser  performance  to  variations  in  uncertain  rate 
coefficients  is  shown  and  rate  coefficient  measurements  critical  to  improved 
modeling  are  indicated.  The  uncertain  rate  coefficients  for  the  more  important 
mechanisms  include: 

1.  Rate  of  deactivation  of  DF(v)  by  D  and  F. 

2.  Overall  pumping  rate  and  distribution  into  vibrational  levels  of 
DF  by  F  +  D£  and  D  +  F2. 
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3.  Rate  coefficient  for  transfer  from  DF(v)  to  C02(00®1)  as  a 
function  of  T  and  v. 

4.  Rate  coefficient  for  collisional  deactivation  of  CO2(00®l)  by  DF 
over  a  wide  temperature  range. 

5.  Rate  coefficient  for  most  of  the  intermolecular  and  intramolecular 
VV  reactions  in  CC>2  involving  levels  (1 1  * 0),  (03*0),  (10°0),  and 
(02*0). 

Additional  measurements  and  theoretical  estimates  of  these  rate  coefficients 
will  improve  the  ability  of  theoretical  laser  models  to  predict  laser 
pe  rformance. 

Many  of  the  recommended  rate  coefficients  are  tentative  and  subject  to 
reevaluation  as  more  data  become  available.  Nevertheless,  the  relationship 
of  pumping,  transfer,  deactivation,  and  lasing  and  the  competition  between 
deactivation  of  OF  and  transfer  to  CO^  has  been  clearly  revealed  by  these 
calculations. 
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